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ABSTRACT 

Simple momentum-transfer theory for the composition dependence 
of ion mobilities and diffusion coefficients in gas mixtures at 
arbitrary field strengths is corrected and extended, and compared 
with a similar theory based on momentum and energy transfer, and 
with results based on direct solution of the Boltzmann equation by 
Kihara’s method. Final equations are recommended for predicting 
composition dependences, given only results on ion mobilities and 
diffusion coefficients in the pure component gases* 
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Simple momentum- transfer theory has recently been used to pre- 
dict the composition dependence of ion mobilities^ and diffusion co- 
2 

efficients in gas mixtures at arbitrary field strengths. The p\ir- 
pose of this comment is threefold: to extend the rcinge of validity 

of the results of Ref. 2 to higher fields by allowing the ion tem- 
perature to be anisotropic and to differ from the gas temperature; 
to correct an error in the expression used for the ion energy that 
affects the final results of Refs. 1 and 2; and to point out the 
relation of these results to a similar theory developed by Milloy 

3 

and Robson, and to recent work based on direct solution of the 

4 

Boltzmann eqxiation. 

By equating the momentum transferred to ions of charge e by an 

electric field, E , to the momentum transferred by the ions to the 

neutrals through collisions,^ and by taking into account the pressure 

forces that arise if the ions are not distributed uniformly in 
2 

space, one obtains the following momentum-balance equation for the 
mean ion velocity ^ v ) , 

5n<v)Ep.N^ ^ ^ ^ ^ » (D 

where 5 is a proportionality constant of order unity, n the 
number density of ions, the reduced mass of an ion-neutral 

pair, Nj the number density of neutral species j , v^^ the 
relative speed of a colliding ion-neutral pair, their 

momentum-transfer or diffusion cross section, and p^ the partial 
pressure tensor for the ions. Pointed brackets indicate an average 
over relative speeds. In a coordinate system in which E is directed 



- 3 - 


along the z-axis, is a diagonal tensor with components 

^xx ~ ^yy ~ ^J_ ^zz ” ^ I ! * that the ions obey the 

ideal gas equation of state, we introduce 


Pjj =nkTjj , Pj. ■ 

where the temperatures parallel and perpendicular to the field are 
defined by 

l^T]! = m [<v^2> - <v^>^J . 

where m is the ion mass. Then V ♦ p^ = kT^ • Vn . 

We focus attention first on the mobility by considering spa- 
tially uniform conditions, for which ( v) = = KE , where 

is the drift velocity and K the mobility. From Eq. (1) 
we can write K in terms of the mobilities K. in the pure com- 
ponents , 


1 . ’<1 


(4) 


where the Xj are mole fractions, and ( means an average in 
pure gas j , as distinguished from an average in the mixture ^ ), 
The evaluation of the ratios {v ,Q. ) / ( v_^.Q.S. proceeds as in 
Ref. 2, and involves three approximations. First, the average of 
the product, ^v^^Q(v^^)^ , is decomposed into the product of the 

averages, )• Second, the average { ^p j ) is re- 

placed by the rms average \ v^. ) ' , which is then evaluated 
from the partitioning of the ion energy among thermal energy, drift 
energy, and random field energy. The third approximation enters 
through the expression for the energy partitioning. 



<v^> + <V^^> = 3kT/n. + <v^>2 [l + <M>^ym]/«) 
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where V. is the neutral species velocity, M. its mass, T 
is the gas temperature, and (m ) is a mean mass of the gas 
mixture, defined as 


<M> 


av 


Eb).M./Eu>. 
j D ] j : 


* 


( I ). 

D 




( 6 ) 

(7) 


The energy partitioning can be calculated by Wannier’s method, 
as was done in Ref. 1, or equivalently by an energy -balance 
equation, as was done in Ref. 3. Unfortunately, there is an error 
in the expression for the energy partitioning in Ref. 1 — in Eqs. 

(27) and (29) the should be replaced by (v^^Q^ ^ to give 

the correct weight factors «. as in Eq. (7). The numerical 
example in Ref. 1 is therefore incorrect in detail because of the 
incorrect used. However, at high fields the relative speed 

is the same as the ion speed, and there is then no mistake in the 

(0 • • 

1 

The theoiy of Milloy and Robson^ for the mobility is essentially 
the same as the foregoing treatment up to this point, but makes 
further approximations in order to express the deviations from Blanc's 
law entirely in terms of the behavior of the ions in the pure neutral 
components, without explicit reference to the mixture itself. ‘Ihus 
the mixture averages aupe e:q>anded in series 


<''rj = <%! j 


d<v . > . 
N r] / 3 


( 8 ) 
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From Eq. (5) we find 


<^/> - j = ''d^ b ^ <H>av''’] - ''dj^ ^ “j'"’ • 


Since this expression appears as a correction term in Eq. (8), it 
need be evaluated only in lowest order, corresponding to 
Qj> ?= <(v^j ^ , which yields 

^ — -T 

^ J J U J ^ V, .(id + M.) 


(m + M.) V,. 
3 dj 


The derivative in Eq. ( 8 ) is obtained by differentiation of Eq. (1) 
written for a single pure component having the same total number 
density as the mixture (N^ = M), and is found to be 



- ^ 2 'm t M. / d In (E/N) 

^ ^dj 3 

Substituting these results back into Eq. (4), we obtain 


, X. - X. d In K. 

F - ^ ^ T ^ ^ ( 1 - A ) , 

j ^ j : d In (E/N) J 


A."^ = (m + M.) K.^ 
3 j 3 


( 12 ) 


( 13 ) 
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which is the result of Milloy and Robson. This expression is re- 
markably similar to one obtained from a direct solution of the 

li 

Boltzmann equation by Kihara’s method. The main difference is that 
the Aj in Eq. (12) is multiplied by a complicated factor that 

happens to be numerically close to unity. Thus Eq, ( 12 ) can be recom- 
mended as an excellent estimate for deviations of the mobility from 
Blanc's law. 


The preceding calculations can obviously also be applied to the 
longitudinal (Dj | ) and transverse (D_l) diffusion coefficients, 
but an additional feature must enter. It is clear from Eq. (1) 
that spatial inhomogeneity affects < v) . since Eq. (5) shows 
that is related to the mean ion speed, spatial inhomo- 
geneity must also affect ^v^^^ . Thus we must allow both 

^ V ) ^ ^rj ) have small contributions from spatial inhomo- 

geneity , but we may consider these as perturbations and linearize 


with respect to them. Denoting the spatially homogeneous case 
with a superscript , we expand the average 
in a series similar to Eq. (8), 


<v .Q(v .)> 
r]^ r] ' 


f<v .2> - 




(14) 


From Eq. (5) we find, on setting {v = v 

' z ^ d ^ 




VdkT 


neE 


3n 

3z » 


(15) 
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the last step following by substitution for from Eq. 

(1). Since this last result is independent of j , we can multiply 
Eq. (14) by over j.; the derivative term is then 

obtained by differentiation of Eq. (1), yielding a result similar 
to Eq. (11) » 


d In ^ 

d {v ° 

' r] ' 


if m 1 d In K 
2v 2 kV-rH-)— I d ln(E/10 

d 


(16) 


Substituting these results back into Eq. (1) and linearizing with 
respect to derivatives of n , we identify the transport coefficients 
and observe the following relations among them: 


eD 


K 


ii 


kT. 



d In K "I 
d ln(E/N) J 


(17) 


eD, 


^ kT, 


(18) 


Equations (17) and (18) have been derived previously by pheno- 

5 6 

menological arguments, * and have been shown to give good quali- 
tative agreement with numerical solutions of the Boltzmann equa- 

6 4 

tion. A more exact treatment according to kinetic theory gives 


eD 


< 1.1 


= k T 


II, J. 


1 + Y. 


d In K 


d in (E/N) 


+ ... 


(19) 


where Y|| 3/2 and ^1/2 , which is valid to order (E/N)^. 

Within this order, the ratio is the same whether 

Yji = 3/2 or 1 and Yj_ = 1/2 or 0 are used. Thus Eqs. (17) - 
(18) are consistent with Eq, (19), even though the Individual 
Y values differ. 

The ion temperatures are obtained from Eq. (3) by an energy- 


partitioning calculation, either by Wannier’s method applied to 
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mixtures , or from a first-order solution to the Boltzmann equa- 
7 « 

tion. In either case, the result for the mixture is 


kX*, = kX t 


5m - (2ni - <M> J) <MA*> 


5m -h 3 


<M*> 


V ^ 


kT_^ = kT + 


5m + 3 ^MA } 


where is given by Eq. (6) and 




S 2 «. M. A. /I (d. , 
j D 3 : j j 


in which A is a ratio of cross sections of order unity • 

To find the composition dependence of the diffusion coefficients 
we substitute Eq. (19), which applies to both mixtures and pure 
gases, back into Eq. (4), 


IlLi 




Tn <v .Q(v .)>° L 


rj rj 


d InK. 

'i 

I ,id in(E/N) 


, ^ d InK 1 

^ll^ld ln(E/N) J 


The important new feature is the appearance of the ion temperatures in 
Eq. (23). The temperatures in the mixture and in the pure components 
are not the same at the same value of E/N because the partitioning 
of the ion energy among drift and random field components depends on 
the nature of the ion-neutral collisions. One immediate consequence 
of this is the recognition that the diffusion coefficients will deviate 


from Blanc's law at high fields even for the Maxwell model of constant 



g 


collision frequency, for which the ratios j 

are all unity and all derivatives of mobility with respect to E/N 
are zero. Numerical calculations of deviations from Blanc's law 
for the Maxwell model have been performed and the results are shown 
in Fig. 1. 

The ratios )j in Eq. (23) can obviously 

be evaluated by the same methods used for the mobility. Substitut- 
ing these results into Eq, (23) and linearizing with respect to the 
derivatives of mobility, in order to be consistent with the trunca- 
tions in Eqs. (8) and (14), we obtain an expression that still con- 
tains a mixture quantity on the right-hand side, namely the deriva- 
tive d InK/d ln(E/N) from the last factor of Eq. (23). In order 
to evaluate this quantity in terms of pure component properties, we 
differentiate Eq. (12), which leaves us with second derivatives of 

Kj with respect to E/N. These second derivatives can be approxi- 

0 

mated by appeal to the following symmetry argument. Since the gas 
is isotropic K must be an even function of E/N; expansion in series 
must therefore involve only even powers of E/N, 


K = K(0) 


jl + Oj^(E/N)^ + ajCE/N)^^ + ... 


(24) 


Differentiating twice and comparing results, we find 


~ — 9 % ^ H tni'F/n) of Order (E/N)** 

d (In E/N)2 ^ 


(25) 


This expression holds for both mixtures and pure gases. Substitu- 
tion of this expression then yields the final result for the com- 
position dependence of the diffusion coefficients, 
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i j: /-3 


I 


X. 

(Dn^l ) 


L 


^'^11,1 ^-i ~ '^HJ. 


M 


2kT| j d In (E/M) * 2Y| | . 


( 26 ) 


The ion temperatures in this expression are to be evaluated from 
Eqs. (20) - (21), but with the average masses calculated with weight 
factors ( 1 ). approximated by making expansions as in Eq. (8), 


X. 

X 


“j (m + H.) K. . 

3 D 1_ 


. X. d In K. 

1 ♦ I ? k7 i in - Ce ThT - ‘i> 

X ^ X * 


(27) 


These results are nearly the same as those obtained by solution of 
the Boltzmann equation. The differences are relatively minor: the 

Aj is missing the factors as in the mobility expression of 

Eq. (22), the values of Y|| and Yj^ somewhat different, and 


the Boltzmann expression does not include the small derivative term 
in the expression for the u. . 

As an excellent estimate for deviations of the diffusion coef- 
ficients from Blanc's law we recommend the following procedure: use 

Eq. (26) with A^ given by Eq. (13), with Yj | and Yj^ given their 
accurate kinetic-theory values, and with the given by just the 

leading term of Eq. (27). 

As a numerical example we consider ions in an equimolar 

mixture of He + Ar at 300® K, using as input data the measurements 
+ . 3 + . 9 

on K in He and K in Ar. Deviations from Blanc's law for K 
and for Djj are shown as a function of field strength in Fig. 2, 
as calculated from the Boltzmann equation and from Eqs. (12) and 
(26). The shapes and magnitudes of these curves have been discussed 

4 

elsewhere ; the iniportant point here is the closeness of the 
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agreement for the two calculation methods. 

In summary, the various theories for the composition dependence 
of ion mobilities and diffusion coefficients in gas mixtures have been 
shown to be consistent. Equations (12) and (26) (modified as 
suggested above) are recommended for predicting deviations from 
Blanc's law as a function of field strength. 
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FIGURE CAPTIONS 


FIG. 1. Ratio of D|j (Blanc) to D||(mix) at high fields for 
an equimolar Maxwell -model binary gas mixture, as contours in a 
^1^^2 have taken m = M 2 • Notice that 

the deviations from Blanc's law are always negative. The corre- 
sponding ratio for mixture mobility is always unity for the Maxwell 
model. 

FIG. 2. Percentage deviations from Blanc's law for mobility, 

100 Ik”^ - iC} / , and for D ,1 as a function of field 

j_ Blanc^ / Blanc M 

strength, as calculated for K'*’ ions in an equimolar He + Ar 
mixture at 300®K. The solid curves are the results from the Boltz- 
mann equation according to Ref. 4, and the dashed curves are the 
present momentum-transfer results of Eqs. (12) and (26). 







